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The resistance of multilayers with interface roughness on a length scale which is large compared 
to the atomic spacing is computed in several cases via the Boltzmann equation. This type of 
roughness is common in magnetic multilayers. When the electronic mean free paths are small 
compared to the layer thicknesses, the current flow is non-uniform, and the resistance decreases 
in the Current-Perpendicular-to-Plane (CPP) configuration and increases in the Current-In-Plane 
(CIP) configuration. For mean free paths much longer than the layer thicknesses, the current flow 
is uniform, and the resistance increases in both the CPP and CIP configurations due to enhanced 
surface scattering. In both the CPP and CIP geometries, the giant magnetoresistance can be 
either enhanced or reduced by the presence of long length scale interface roughness depending on 
the parameters. Finally, the changes in the CPP and CIP resistivities due to increasing interface 
roughness are estimated using experimentally determined parameters. 

PACS numbers: 75.70.Pa, 73.40.-c 



I. INTRODUCTION 

The study of metallic multilayers has been a very ac- 
tive area of research in recent years. Of particular in- 
terest are alternating layers of ferromagnetic and param- 
agnetic metals called magnetic multilayers. A relatively 
small magnetic field aligns the magnetic moments in the 
different ferromagnetic layers, leading to a large mag- 
netoresistancej-,|Wliich is called the giant magnetoresis- 
tance (GMR) Ma The GMR has technological applica- 
tions in magnetic read heads, magnetic sensors, and mag- 
netic memory devices. Consequently, it has been studied 
extensively, including such effects as bulk and interface 
scattering, magnetic and non-magnetic scattering, Fermi 
wave vecto^-mismatch and magnetic coupling between 
the layers Boflla 

One aspect of metallic multilayers which has been less 
studied theoretically is the effect of long length scale fluc- 
tuations of the layer thicknesses and heights. In an ideal 
multilayer the interfaces between the layers are perfect 
planes. The thickness of each layer and the height of 
each layer above the substrate would be constant. Obvi- 
ously, this is not the case in any real system. In addition 
to interdiffusion and other atomic scale disorder at the 
interfaces, the actual thicknesses and/or heights of the 
layers can vary on a length scale which is large on the 
atomic scale. Indeed, these long length scale fluctuations 
seem to be the rule rather than the exception. 

Since the Fermi wave vectors for these metals are of 
order the atomic spacing, interdiffusion and atomic scale 
disorder are strong sources of scattering. Fluctuations on 
a length scale of 10 or more atomic spacings should not 
appreciably effect the surface scattering. Nonetheless, 
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long length scale disorder can be important. It has been, 
demonstrated both experimentallyu and theoretically^ 
that nonplanar interfaces can create new kinds of mag- 
netic coupling between the layers. This tends to reduce 
the GMR because the fraction of a sample that is anti- 
ferromagnetically aligned at zero applied field is reduced. 
Consequently, many of the experiments that study the 
role of interfacial roughness rescale the magnetoresis- 
tance by the fraction of the sample which is antiferro- 
magnetically coupled. This fraction can be determined 
experimentally from magnetization measurements. 

There are two common geometries used in studying 
the giant magnetoresistance. The current flows paral- 
lel to the layers in the Current-In-Plane (CIP) geometry 
and normal to the layers in the Current-Perpendicular- 
to-Plane (CPP) geometry. The CIP geometry has been 
more widely studied than the CPP geometry because 
measurements in the CIP configuration are easier to 
achieve experimentally. Simultaneous measurements of 
the CIP magnetoresistance and roughness using low an- 
gle X-ray scattering have been made on samples with 
roughness that has been systematically varied by chang- 
ing growth parameters and annealing. One set of ex- 
periments on Fe/Cr multilayers finds that the change in 
the resistivity between low and high field, Ap, rescaled 
by the antife r.mmagnetic fraction, AFF, increases with 
roughness .oE£rEil Schad et al. also simultaneously mea- 
sured the surface rotjghness and the magnetoresistance 
of Fe/Cr multilayers. E3E3lj In one set of experiments 
they find that Ap/AFF decreases with the increasing 
fluctuations in the-iasier heights while the saturation re- 
sistivity decreases. EiEj In another set of experiments on 
Fe/Cr multilayers which are dominated by surface scat- 
tering they find that the CIP magnetoresistance increases 
with interface fluctuations, which were determined by the 
ratio of the vertical, roughness amplitude to the lateral 
correlation length.EJ 

The effect of long length scale roughness on the CIP 
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magnetoresistance has been studied theoretically by Bar- 
nas and Bruynseraeddi3 and by Levy et al£3 Barnas 
and Bruynseraede studied the scattering between quan- 
tum states in different layers allowing for uncorrelated 
variations in the layer thicknesses. They find that the 
magnetoresistance can either increase of decrease with 
roughness depending on the parameters in the prob- 
lem. Indeed there are even cases where the magnetoresis- 
tance decreases and then increases with roughness. Levy 
et al. studied the magnetoresistance of multilayers de- 
posited on grooved substrates J13't3 which produce rough- 
ness which is correlated between the layers. Using a gen- 
eral linear response approach they find that the rough- 
ness mixes the CIP and CPP geometries, leading to what 
is called the CAP or Current-at-an-Angle-to-Plane con- 
figuration. Since the CIP magnetoresistance is typically 
smaller than the CPP magnetoresistance, this would tend 
to increase the magnetoresistance with roughness. 

The difficulty in making CPP measurements is due to 
the large surface area relative to the thickness of the sam- 
ples. This results in a CPP resistance which is too small 
to measure using conventional techniques and also makes 
the measurements susceptible to inhomogeneous current 
paths. These difficulties hajie been overcome by mak- 
ing small area pillar samples Jij by using superconducting 
leadsjH3 and by combinations of these two techniques. Ell 
Experiments on Co/Ag multilayers find that the CPP 
magnetoresistance decreases with interface disorder while 
the CIP magnetoresistance increases.Ea Cyrille et al. mea- 
sured the CPP magnetoresistance of Fe/Cr multilayers 
and quantified the roughness in their samples using both 
low angle X-ray scattering and trans m i ss ion electron mi- 
crographs of cross-sectional samples. EilO They find that 
Ap/AFF increases with roughness proportionally to the 
RMS deviation in the layer heights. Still more recent ex- 
periments by Zambano et al. see no change in the CPP 
magnetoresistance of Fe/Cr multilayers with increasing 
roughness .cil 

The effect of long length scale interface fluctuations is 
therefore unclear, with some work pointing to an increase 
in the magnetoresistance, some a decrease, and some no 
change at all. Some of the differences between the ex- 
periments is due to changing a growth parameter like 
the sputtering pressure probably changes more than just 
the interfacial roughness. A theoretical calculation can 
vary just the interface fluctuations and hence hope to iso- 
late the effect of long length scale interface fluctuations 
on the magnetoresistance. While the earlier theoretical 
work is consistent with the CIP magnetoresistance exper- 
iments, the CPP magnetoresistance experiments remain 
unexplaiaed. In particular the extensive work of Cyrille 
et alr3T3 which shows an increase in the magnetoresis- 
tance with interface roughness is not possible to explain 
as a mixing of the CPP and CIP geometries since the CIP 
magnetoresistance is usually lower than the CPP one, as 
is the case in their measurements. 

In this paper we examine the effects of long length scale 
fluctuations in the layer thicknesses and heights using the 



Boltzmann equation. A semiclassical approach like the 
Boltzmann equation is a good choice for this problem 
because the length scales involved are large compared to 
the atomic spacing. Various versions of the Boltzmann 
equatkm have also been used extensively in modeling the 
GMR.EE3t3'E!l The Boltzmann equation we use here has 
a simple, but current conserving form for the scattering 
term, and we can solve it essentially exactly numerically 
in the limits of long and short mean free paths. Here, a 
long mean free path is much larger than the layer thick- 
nesses, and a short mean free path is much smaller than 
the layer thicknesses. We consider both the CPP and 
CIP geometries and find that the GMR can either in- 
crease or decrease with interface roughness depending on 
the parameters in the problem. Explicit predictions are 
made for when the GMR increases and when it decreases. 

The remainder of the paper is organized as follows. In 
Sect. II we provide a formal solution to the Boltzmann 
equation which is valid for arbitrary mean free paths. 
Next, in Sect. Ill we apply this formal solution to repro- 
duce the well known results for the CIP and CPP geome- 
tries when the interfaces are flat. In Sect. IV we consider 
sinusoidal interfaces, which model the long length scale 
interface fluctuations. Both the long mean free path case 
(Sect. IV A) and the short mean free path case (Sect. IV 
B) are considered. These two cases are compared numer- 
ically in Sect. IV C. In Sect. V we use the results from 
Sect. IV to determine the effect of long length scale inter- 
face disorder on the giant magnetoresistance in the long 
mean free path case (Sect. V A) and the short mean free 
path case (Sect. V B). In Sect. V C the changes in the 
CPP and CIP resistivities due to increasing long length 
scale roughness are estimated in the long and short mean 
free path cases using experimentally determined param- 
eters. All the results are summarized in the conclusion, 
Sect. VI. 



II. FORMALISM 

The model system we study is shown in Fig. 0. Mul- 
tilayers of thickness Ay are separated by sinusoidal in- 
terfaces with amplitude A and period £. Within each 
layer the relaxation time, Tj, is constant. For numeri- 
cal considerations, our calculations are performed in two 
dimensions, namely the multilayers are strips in a two- 
dimensional plane; however, we nonetheless expect the 
results to remain qualitatively the same when general- 
ized to three dimensions. Moreover, in both the long 
and short mean free path limits we will present analytic 
expressions for our results, which have direct generaliza- 
tions to the three-dimensional case. 

The Boltzmann equation we use represents elastic s- 
wave scattering within a current-conserving right-hand 
side, 



v ■ V r f - eE • V p / = - 
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FIG. 1: Schematic of a four-layer repeat unit with interfacial 
roughness. The interfaces are modeled as sine waves of ampli- 
tude A and period £. Layer i has a thickness Ayi and uniform 
relaxation time n, and the total thickness of the repeat unit 
is Ay. We calculate the current density and conductivity due 
to roughness for currents flowing in the y-direction (CPP con- 
figuration) and ^-direction (CIP configuration) for mean free 
paths which are long and short compared to Ay. 



where / = /(r, p) is the distribution function, / = 
/(r, |p|) is the spherical average of / in momentum space, 
and t — r(r) is the relaxation time. To get the linear re- 
sponse conductivity the distribution function is expanded 
to first order in the applied electric field, / = f eq + Sf, 
where f eq is the equilibrium Fermi-Dirac distribution 
function and 6f is proportional to the applied electric 
field, E. The linear response Boltzmann equation is 



v ■ W r Sf - eE • v 



\ de 



^ Sf-Sf 



(2) 



This equation may be further simplified using the fact 
that the temperature is far below the Fermi tempera- 
ture in these metals. Consequently, the energy derivative 
of the equilibrium Fermi function in Eq. (H) is approx- 
imately a delta function which pins the energy to the 
Fermi energy. One can define another distribution func- 
tion on the Fermi surface, g, 



Sf 



Of 



eq 



de 



so that the Boltzmann equation becomes 



v • W r g + eE • v 



(3) 



(4) 



This first order differential equation can be integrated 
to find g and hence /. As a first step the deviation of g 
from its spherical average ~g is defined as Sg, 



Sg = g- 9- 



(5) 



Substituting 5g into Eq. (0) , the Boltzmann equation on 
the Fermi surface is now 



W r 5g 



Sg 



-eE v v V,..g. 



(G) 



Both the gradient of the electrical potential, E = — V0, 
and the gradient of g appear on the right hand side of 
Eq. (||). Consequently, it is useful to define the electro- 
chemical potential, V, 



V = <j>- 
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so that there is only a single gradient on the right-hand 
side, 



v • V r <5<7 



Sg 



= ev • W. 



(8) 



Equation (|8|) can be solved by integrating along paths 
in phase space. Letting (r(s),p(s)) be a trajectory in 
phase space which satisfies the equation of motion, f = v 
and p = 0, the distribution function along this path is 



g(s) = 6g{r(s),p(s)). 



(9) 



Substituting Eq. (g) into Eq. (g), one obtains a first order 
ordinary differential equation, 



ds 



g _ dV 
t ds 



The general solution of this equation is 

«/ i 



g( s f) 



exp 

rs f 



exp 



t(s) 
f 



dsj g(si) 

) 6 ^ 



(10) 



(11) 



where Sj and s / represent the initial and final coordinates 
in phase space, respectively. Note that with these sim- 
ple equations of motion for this Boltzmann equation, the 
trajectories are lines, p(s) = p; = P/ and 



r ( s ) = r / - v ( s / - s )- 



(12) 



When the starting point of integration is taken to infinity, 
Si — > — oo, Eq. (fill) simplifies to 



exp 



1 



V) 



-ds' 



e——ds. 

ds 



(13) 



Equation (|13|) is the starting point for the calculations in 
this paper. Solving it is not always straightforward be- 
cause the electro-chemical potential is not known a priori. 
In some cases one can deduce V from general principles. 
In other cases Eq. (|l^) must be solved self-consistently 
as described below. 

Once one has a solution to g = Sg(r(s),p(s)), the cur- 
rent density, j, can be obtained directly from 

r J/1 

Ja (r) = -eN(E F ) / ^ Va 5g(r,8 p ) (14) 
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because g has no angular dependence. In this equation 
the density of states at the Fermi surface is N(Ep) = 
m/(irh 2 ). For the nonuniform samples considered in this 
paper, the current density and the electric field will vary 
with position. An effective conductivity, a, for the entire 
sample is defined using the spatial average of the current 
density, (j), and the electric field, (E), 

{3a)=CTaA E p)- (15) 

This is the conductivity one would obtain by measur- 
ing the net current and voltage drop for a large sample 
and multiplying by the usual factors of length and cross- 
sectional area. In Eq. ( |l5|) the conductivity is defined in 
terms of the electric field instead of the gradient of the 
electro-chemical potential. However, because of the pe- 
riodicity in the multilayers, ~g is periodic. Thus (Vg) = 
because the spatial average of the derivative of a periodic 
function is zero. Consequently, from Eq. (FtI) the average 
electric field, (E), is the same as the average gradient in 
the electro-chemical potential, (— VV). 



III. FLAT INTERFACES 

In this section we apply the formalism developed in the 
previous section to solve the case of flat interfaces (A = 
in Fig. |l|) . This will provide the basis for the calculations 
with wavy interfaces in the next section, since in this 
case we are able to solve the Boltzmann equation exactly 
for arbitrary mean free paths. It will also allow us to 
demonstrate that this formalism reproduces the conven- 
tional results in the limit of weak spin-flip scattering. 



A. CPP geometry 

In the Current-Perpendicular-to-Plane (CPP) geom- 
etry the current flows in the y-direction of Figure [l]. 
Current conservation and translational symmetry in the 
x-direction imply that the current density is uniform 
throughout the sample. Translational symmetry also im- 
plies that the electro-chemical potential, V, only depends 
on the y variable. The actual functional form of V(y), 
however, is not known a priori. 

In order to determine V(y) or equivalently dV /dy, the 
condition that the spherical average of Sg is zero (see Eq. 
(||)) was used in conjunction with Eq. (|l3|). This leads 
to an integral equation for dV /dy. Because of the peri- 
odicity of the multilayers and hence dV/dy, the infinite 
integral can be converted to a finite integral. Discretiz- 
ing dV/dy then leads to a linear equation, which is easily 
solved numerically. 

The results of this numerical calculation show that for 
all mean free paths the derivative of the electro-chemical 
potential in the i— layer is proportional to the inverse 
of t,. Substituting this result into Eq. dl3|), it follows 
that the distribution function is independent of position 



and proportional to y • p, where p is the direction of the 
momentum. The current is uniform as expected. Cur- 
rent uniformity and the fact that dV/dy within a layer 
is proportional to 1/t; imply that the resistance is the 
same as one would obtain by adding classical resistors in 
series. The series resistor model for the CPP geometry 
is commonly used to analyze experiments in the limit of 
weak spin-flip scattering.Ea 



B. CIP Geometry 

In the Current-in-Plane (CIP) geometry the current 
flows in the x-direction of Figure [j]. In this case the elec- 
tric field must be constant because of the symmetry of the 
problem and V x E = 0. Furthermore, if we make the 
ansatz that ~g = 0, then Eq. (|l3| ) can be integrated analyt- 
ically. The spherical average of the resulting distribution 
function, dg, is indeed zero, which is consistent with the 
ansatz that g = 0. Thus, in this case we also have a 
complete analytic solution for the distribution function. 

Figure || shows the CIP current density in a two-layer 
repeat unit multilayer with Ayi = 3 and Aj/2 = 7. Each 
curve corresponds to different mean free paths, ranging 
from much less than to much greater than Ay. For com- 
parison, the curves are normalized by the maximum of ] x . 
As the figure illustrates, when the mean free paths are 
much less than Ay, j x changes rapidly at the interface. 
Moreover, in this limit the ratio of the current density 
in layer 2 to that in layer 1 equals the ratio of the cor- 
responding mean free paths, which is what one would 
expect for classical macroscopic resistors in parallel. As 
the mean free paths increase, the current density smooths 
out, becoming roughly uniform across the repeat unit as 
the mean free paths become much greater than the layer 
thicknesses. 



C. Summary 

In the long mean free path limit the current density 
and the distribution functions in both geometries are in- 
dependent of position. As seen in Eq. (|l3|), there is an 
average over a mean free path which goes into computing 
the distribution function. If the mean free path is large 
enough, then the average effectively smears out the vari- 
ations in the sample, producing a homogeneous current 
and distribution function. 

In the short mean free path limit the current density 
is the same as one would obtain classically from adding 
resistors in series (CPP geometry) or in parallel (CIP 
geometry). In this case the integral used to calculate the 
distribution function, Eq. (|l3]), is short ranged and in 
particular much smaller than the layer thicknesses. Each 
layer behaves like a macroscopic piece of metal. 
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FIG. 2: CIP current density for different mean free paths in 
a two-layer repeat unit with Aj/i = 3 and A3/2 = 7. The mean 
free path in layer i is U — VFTi. For comparison, the current 
density for each curve is normalized by its maximum value, 
max(ja;). In the short mean free path limit (£1,^2 -C Ay) the 
current density changes sharply in proportion to £i/t2 across 
the boundary. Conversely, in the long mean free path limit 
(£i,£2 Ay) the current density is roughly uniform. 



IV. CURVED INTERFACES 

For the case of the wavy interfaces shown in Fig. [j] 
it is not possible to obtain simple analytical expressions 
for the distribution function like the ones for the flat in- 
terface cases of the previous section. Moreover, direct 
numerical solution of the Boltzmann equation requires 
the solution of a three-dimensional problem (two space 
coordinates and one angular variable along the Fermi sur- 
face). However, we saw in the previous section that there 
are two natural limiting cases: (i) the long mean free path 
case, where the current density and distribution func- 
tion are uniform, and (ii) the short mean free path case, 
where the current density and distribution function are 
determined by the local gradient in the electro-chemical 
potential. In this section we solve for the distribution 
function, current density, and conductivity in these two 
limiting cases. 



A. Long mean free path limit 



In the case where the mean free path is long compared 
to the layer thicknesses the integral used in computing 
the distribution function (Eq. |l3|) averages the gradient 
in the electro-chemical potential over a large region of the 
sample. If the mean free path is long enough, then the 
gradient in the electro-chemical potential may be approx- 
imated by its average value. The distribution function 



then becomes 

g( s f) 



-ev(E)g'( S/ ) 



exp 



>t 1 



r(s') 



ds' I ds, 



(16) 
(17) 



where the angular brackets denote a spatial average. As 
discussed in Sect. II, the average of the gradient in the 
electro-chemical potential, (— VV), is equal to the aver- 
age electric field, (E) , for a periodic system. 

In the previous section, where we considered flat in- 
terfaces, we did not include interface scattering in the 
interest of simplifying the calculation. The scattering 
rates changed as electrons went from one material to an- 
other, but there was no additional scattering due to the 
interfaces. Here we include interface scattering since it 
is crucial for comparing to experiments. We model a 
simplified surface scattering as infinitesimal layers of a 
higher-resistivity material. Surface scattering is included 
on the right-hand side of the Boltzmann equation of Eq. 
([!]) in addition to the bulk scattering term, 1/rn,, 



1 



1 



r(r) r 6 (r) 



r]T y"dM 2 (r-RK0), 



(18) 



where the integration runs along the i— interface, Ri(l) 
is the position of the i— boundary, and T is a parameter 
characterizing the surface scattering rate. Letting £ be 
the angle between the velocity, v, and the unit vector 
normal to the boundary with the convention that < 
C — 7r /2, the integral in Eq. ( |l7| ) from slightly below 
(s = t-) to slightly above (s = t + ) the interface is 



1 



t(s) 



ds 



vp cos(C) ' 



(19) 



Physically, this means that the probability of an electron 
being scattered at an interface is lowest when its trajec- 
tory is perpendicular to the boundary (( = 0) and larger 
when it crosses the boundary at an angle. 

To test the approximation of Ens. (|l^) and (|l7|) we 
now evaluate the integral of Eq. ( |l7| ) numerically. The 
resulting function, g' , should be independent of both 
angle and position. Figure I shows g'{6) versus 9 for 
a multilayer composed of two-layer repeat units with 
£ = 10, Ayi = Ay 2 = 2.5, l x = 1000, i 2 = 2000, and 
Wi?/r = 500. The dashed and solid lines represent this 
sample at A = and A — 1, respectively. The statistical 
noise present in the latter curve near 9 = and 9 = tt 
corresponds to electrons whose trajectories run nearly 
parallel to the x axis. Although these electrons intersect 
significantly fewer interfaces than those that travel along 
the y axis, the interfaces can nevertheless be much more 
effective in scattering them if they intersect a boundary 
nearly tangentially within a few mean free paths. There- 
fore g' is highly sensitive to small changes in 9 near 9 = 
and 9 = tt. Despite the magnitude and frequency of 
the fluctuations in this region, the noise averages to the 
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FIG. 3: The function g'(8) used to compute the distribution 
function in the long mean free path limit. The fluctuations 
present near 6 = and 8 = n are due to paths which do not 
intersect many interfaces. Once these rapid oscillations are 
averaged locally, however, this function is roughly indepen- 
dent of angle and position as expected in the long mean free 
path limit. The decrease in g' as one goes from flat interfaces 
(A = 0) to curved interfaces (A = 1) leads to a decrease in the 
interface conductivity with roughness in both the CPP and 
CIP cases. The data shown correspond to a two-layer repeat 
unit with £ = 10, Ayi = Ay 2 = 2.5, t x = 1000, l 2 = 2000, 
and v F /V = 500. 



smooth value of g' seen in the vicinity of 9 = ±7r/2. This 
kind of averaging is performed to c omp ute the current, 
Eq. (|l4|), and the conductivity, Eq. (15). Thus, once the 
small angle fluctuations are averaged the function g' is 
independent of angle. We have also verified that it is 
independent of position for a given sample. 

From Fig. |j| it is apparent that g' decreases with in- 
creasing interface roughness (A), resulting in a decreas- 
ing conductivity in both the CIP and CPP geometries. 
Performing the same calculation except without surface 
scattering (r — 0), one finds that the two values for g' 
are the same. In other words, with only bulk scattering 
in our model interface roughness does not effect the con- 
ductivity. In the long mean free path limit, it is therefore 
necessary to have interface scattering in order to have a 
change in the conductivity with roughness. 

In order to examine the additional contribution to the 
resistance due to interface scattering, we define an in- 
terface conductivity, a*, which is extracted as follows. 
For a given geometry the conductivities with and with- 
out interface scattering are computed, cr^o an d 0T=o, 
respectively. Treating the bulk and interface resistances 
as resistors in series, the interface conductivity is given 
by 



(. 



M 



Ay \ ar^o <Jr=o 



(20) 



is the thickness of a repeat unit. In the example above 
the number of layers is M — 2 and the length of a repeat 
unit is Ayi + A 7/2 = 5. Using a resistors in series model 
is natural in the CPP geometry; however, it may be less 
clear that this is a good model for the CIP geometry. 
Here, because g' is independent of angle, the changes in 
the CPP and CIP conductivities are equal. Thus, what 
works in one geometry will also work in the other. 

For a given set of bulk and surface scattering rates 
define 8a* as the difference between the interface con- 
ductivity for A 7^ and A = 0, 



6a* 



a*(A)-a*(0). 



(21) 



Figures [| (a) and (b) contain plots of \5a*/a*\ versus 
(A/£) 2 for the CPP and CIP geometries, respectively. 
From Eq. ( |l6|) one can see that the distribution function 
Sg — g in the CIP case, where the average electric field is 
in the rr-direction, is largest near 9 — and 9 = it, while 
in the CPP case, 8g is largest near 9 = ±7r/2. The larger 
scatter in the data of Fig. [| for the CIP case than the 
CPP case reflects the larger scatter in g' near 9 = and 

9 = 77. 

The curves in these figures represent either various 
positions within a particular sample or fixed positions 
within samples that have different mean free paths or 
surface scattering rates. When plotted in this manner, 
all of the data fall onto the same line. Thus, not only is 
da* /a* proportional to (A/£) 2 , but the proportionality 
constant is independent of the model parameters. 

One difference between a curved interface and a flat 
interface is that the curved one is longer. Because of sur- 
face scattering, a curved interface will have more scatter- 
ing and hence a larger resistance. The amount of extra 
scattering provided by a curved interface should be pro- 
portional to the additional length in the boundary. Let 
C be the length of the interface and 6C be the change 
in the interface length from the flat case (A = 0) to the 
wavy case (A ^ 0). The percent increase in the inter- 
face length, 6C/C, is plotted as the x's in Fig. ^. To 
a good approximation it is evident from Fig. |] that the 
percent decrease in the interface conductivity is equal to 
the percent increase in the length of the interface. For 
our sinusoidal boundaries the percentage change in the 
interface length is tt 2 (A/^) 2 for small A, so within our 
model we have 



6a*_ 
~a* 



6C 



(22) 



where M is the number of layers in a repeat unit and Ay 



B. Short mean free path limit 

In the short mean free path limit the integral used to 
compute the distribution function, Eq. ([l3]) , samples the 
gradient in the electro-chemical potential over a short dis- 
tance. Thus, we may approximate the electro-chemical 
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FIG. 4: Fractional change in the (a) CPP and (b) CIP interface conductivities with roughness in the long mean free path 
limit. Also shown is the corresponding change in interface length (x) with roughness. In both cases the fractional change in 
the interface conductivity, \5a* /a*\, is very close to the fractional increase in the length of the interface, SC/ C. In our model, 
which assumes sinusoidal interfaces, both of these quantities are proportional to (A/!;)' 2 . The (*), (o), and (o) data points 
represent different coordinates within a sample that has Ayr = A1/2 = 2.5, t\ = 1000, £2 = 2000, and vf/T = 500. The (★) 
and (<) points represent the same coordinate as the (o) points, but correspond to samples with £\ = £2 — 1000 and twice the 
surface scattering rate, respectively. The larger scatter in the data for the CIP case is due to the fluctuations in g' in Fig. 0. 



potential as constant within the integral, and the distri- 
bution function becomes 

Sg{T,p)merv F p-W(T). (23) 

In this approximation, the current density and conduc- 
tivity in the i-th layer are given by 

j(r) = <r(r)(-W(r)), (24) 



Since the conductivity is constant within each layer, the 
current conservation condition, V • j = 0, within a given 
layer implies that — V 2 Vi — 0, where Vi is the electro- 
chemical potential in layer i. The problem of finding the 
current density and effective conductivity due to rough- 
ness thus reduces to solving Laplace's equation within 
each layer subject to the boundary conditions that the 
current is conserved and that the electro-chemical poten- 
tial is continuous. As noted earlier, the average electric 
field, (E), is the same as the average of the gradient in 
the electro-chemical potential, —W, because of the pe- 
riodicity of the multilayers. 

Whereas in the long mean free path limit we had to 
include interface scattering or else there would not have 
been any effect, in the short mean free path limit there 
is already an effect without surface scattering. From cal- 
culations which include both bulk and surface scattering, 
we have found that the presence of interface scattering 
does not qualitatively alter the results for the conduc- 
tivity. Thus, in the results presented below we do not 
include interface scattering. 




5 10 5 10 

x x 

FIG. 5: CPP current density vector field and current field 
lines in the short mean free path limit. Due to roughness, 
current tends to traverse a greater distance in high conduc- 
tivity than low conductivity layers, leading to an increase in 
the effective CPP conductivity with roughness. 



Figure || shows the CPP current density vector field 
and current field lines in a repeat unit with A = 1, Ay = 
£ = 10, Aj/j = 2.5, and Oi — io\. As the figure illustrates, 
roughness causes current to flow nonlinearly through the 
sample in such a way that it traverses a greater distance 
through the high conductivity layers than through the 
low conductivity layers. The effect of roughness is thus 
to increase the effective CPP conductivity, a, of Eq. JTe|). 

As in the long mean free path limit, the percentage 
increase in a relative to the flat interface case is propor- 
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^/Ay 

FIG. 6: The proportionality constant acpp between 8a /a 
and (A/ 1;) 2 in the CPP geometry. This constant acpp in- 
creases with the ratio of the layer conductivities and (,/Ay, 
saturating for o\ S> 02 and when £ becomes several times Ay. 

tional to (A/£) 2 , 



Here, however, the proportionality constant, which we 
denote as acpp for the CPP case, depends on both the 
geometry and the layer conductivities. Figure [6] shows 
cicpp as a function of £/ 'Ay for several values of <Ji/<J2 
for a sample composed of two-layer repeat units with 
Ayi = Ay/2. As the figure illustrates, acpp increases 
with the ratio of the layer conductivities and saturates 
for G\ 3> <72. It also increases with £/ Ay, saturating as £ 
becomes several times Ay. 

To determine the current and conductivity in the CIP 
geometry, the same Laplace's equation is solved with the 
net potential drop in the x as opposed to the y direc- 
tion. Figure || shows the CIP current density vector field 
and current field lines in a four-layer repeat unit with 
A = 1, Ay = £ = 10, Ayi — 2.5, and Oi = lo\. Due to 
roughness, electrons in high conductivity layers impinge 
on lower conductivity layers near the interfaces. This in 
effect reduces the short circuit effect of the high conduc- 
tivity layers and consequently decreases the effective CIP 
conductivity. 

The percent change in the effective conductivity rel- 
ative to the flat case (A = 0) is again proportional to 
(A/£) 2 . The proportionality constant, acip, is simply 
related to acpp in the two-layer case: acip approxi- 
mately equals the negative of acpp, where the latter 
is calculated with the layer conductivities interchanged 
{acip{(7u<rz) ~ -atcpp(cr2,(7i))- Equivalently, a C ip 
is roughly equal to the negative of acpp when the lat- 
ter is computed for a sample in which the thicknesses 
of layers 1 and 2 are interchanged (acip (Ayi, Ajfe) ~ 
—acpp(Ay 2 , Ayi)). This result is illustrated in Fig. 




10 20 30 40 

FIG. 7: Comparison of a for the CPP and CIP geometries. 
The numeric CPP and CIP curves shown are roughly mir- 
ror images of one another, illustrating that ctcip(ai, a?) w 
-Qcpp(ff2,ffi). Depending on the sample geometry, the an- 
alytical approximations for a in Eqs. (|2ij) and ( |32[ ) are excel- 
lent in some cases and only roughly correct in others. The 
data shown are for a two layer repeat unit with £ = 10 and 
Ay = 5.2. 
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x x 

FIG. 8: CIP current density vector field and current field 
lines in the short mean free path limit. As a result of rough- 
ness, current near an interface does not remain within a single 
layer across a period but traverses both high and low conduc- 
tivity regions. This reduces the short circuit effect of the high 
conductivity layers and leads to a decrease in the conductivity 
of the multilayer. 

0, where the numeric oicpp versus <J\/(J2 curves for 
Ayi/ Ay = 0.5 and 0.1 are mirror images of the numeric 
o-Cip versus a\/oi curves for Ayi/ 'Ay — 0.5 and 0.9. 

Approximate analytic expressions for both olcip and 
&cp p can be derived-using the approach of Levy et al. for 
grooved multilayerslij. These multilayers are similar to 
those shown in Fig. [l] except that the interfaces are piece- 
wise linear. The angle the interfaces make with the x-axis 
is defined to be 9, and it is assumed that the layers are 
far apart relative to the amplitude of the interface fluctu- 
ations, which we have called A in the sinusoidal case. In 
the CIP configuration they approximate the electric field 



9 



to be uniform in the x-direction. Letting acip and acpp 
be the CIP and CPP conductivities when the interfaces 
are fiat, it follows that the spatial average of the current 
density is 

Ox) = {? C ip cos 2 9 + a C pp sin 2 6) E x . (27) 

Consequently, the effective conductivity is 

a = a cip cos 2 9 + a cpp sin 2 9. (28) 

Although we consider sinusoidal as opposed to piecewise 
linear interfaces, we can use an approximate effective 9 
of 2A/(£/2) = 4A/£. Expanding Eq. (||) for small 9 and 
using the definition in Eq. (^), it follows that a in the 
CIP geometry is 



f &CPP . 

acip s* 16 1 



(29) 



Figure [7] shows acip versus 01/02 for different values 
of Ayi/Ay as computed analytically via Eq. ( p9| ) and nu- 
merically in the short mean free path limit. The data cor- 
respond to a sample composed of two-layer repeat units 
with £ = 10 and Ay — 5.2. As the figure illustrates, the 
agreement between the analytic and numeric acip is ex- 
cellent for some parameters and worse for others. In the 
case of two-layer repeat units, we find the agreement to 
be particularly good when the geometric parameters sat- 
isfy £/Ayi s» 4 for a\ > 02 and £/Aj/2 ~ 4 for <j\ < 02. 
This relation is satisfied by the CIP curves correspond- 
ing to Ayi/Ay = 0.5 in Fig. |^. In any case, the analytic 
expression does provide an estimate for ucip- 

In the CPP case if one starts from Eq. ( p8| ) to obtain 
&CPPi one obtains a poor fit to the numerical solution. In 
the CIP geometry for the case of flat interfaces, symmetry 
required that the electric field be uniform. On the other 
hand, in the CPP geometry for the case of flat interfaces, 
symmetry required that the current density be uniform. 
Thus, in the CPP geometry we start from the assumption 
that the current density is uniform in the y-direction. 
The average electric field in the y-direction is then 



(E y ) = (pcpp cos 2 9 + pc ip sin 2 9) j y , 



(30) 



where pcpp and pcip are the CPP and CIP resistivities 
corresponding to flat interfaces. The effective resistivity 
is therefore 



P = Pcpp cos 9 + pcip sin 9 



(31) 



Note-that this equation is not the same as Eq. ( f28| ) from 
Rcf.EJ, but rather its generalization to the case of the 
CPP geometry. As will be shown below, it provides a 
much better to fit to our numerical data on the CPP 
resistivity. Equation ( |3~l"| ) and the assumption that 9 m 
4A/£ imply that the CPP proportionality constant is the 
negative of the CIP proportionality constant, 



independent of the sample geometry. Note, however, that 
numerically Eq. (|32| ) holds for two-layer repeat units only 
when Ayi = Aj/2- 

The analytic expression for acpp is compared to the 
numerical solution in Fig. |?]. As for the CIP case, the 
analytic estimates are excellent for certain geometries but 
not as good for others. We find the agreement in the 
case of two-layer repeat units to be particularly good 
when the geometric parameters satisfy £/Ay2 ~ 3.7 for 
<j\ > (T2 and £/Ayi w 3.7 for a\ < 02 . The CPP curves 
corresponding to Ayi/Ay = 0.5 in Fig. [7] satisfy this 
relation. 



C. Comparison of the long and short mean free 
path limits 

The physics of the change in the effective conductivity 
due to long length scale interface roughness is different 
in the long and short mean free path limits. The effect 
in the long mean free path limit is due entirely to en- 
hanced interface scattering with roughness, and no effect 
is observed when surface scattering is ignored. In the 
short mean free path limit, interface scattering plays a 
less dominant role, as the effective conductivity changes 
with roughness even when only bulk scattering is present. 
Here the effect in the CPP geometry results from rough- 
ness providing a less-resistive, non-linear path of current 
flow, while the effect in the CIP geometry is due to a re- 
duction of the short circuit effect of the high conductivity 
layers. 

In the long mean free path limit, the decrease in the 
interface conductivity due to roughness approximately 
equals the corresponding increase in the interface length 
in both the CPP and CIP geometries. Thus, in our model 
we have 



5a*_ 
~a* 



(33) 



In the short mean free path limit, the percent change in 
the effective conductivity due to roughness is also pro- 
portional to (A/£) 2 , 



6a 
a 



(34) 



acpp 



-acip, 



(32) 



Here, however, the proportionality constant a depends 
on the layer conductivities and the geometry of the sam- 
ple. For a multilayer composed of two-layer repeat units 
with £ = 10 and Ay\ = Ay2 = 2, the coefficient \a\ is 
approximately 2 for CT1/CT2 = 2, 12 for CT1/CT2 = 10), and 
18 for a\ ^S> (J2, indicating a range of order 10. 

Within each limit, we make predictions concerning dis- 
tinct physical quantities - interface conductivity versus 
effective conductivity. To compare the magnitude of the 
percentage change in the effective conductivity in both 
limits, we write 5a j a in the long mean free path limit in 
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terms of 5a* /a* , where 



1 M 1 



, (35) 
^a r =o Ay a* 1 

Using Eq. (|35|) and the fact that o"r=o is independent of 
A, we find that 

5a 5a* 



Aya* 
Mor=c 



(36) 



Since the second term in the parenthesis is greater than 
or equal to 0, 5a j a must be less than or equal to 5a* /a* , 
where equality holds in the limit of zero bulk resistance. 
Additionally, 5a /a goes to zero in the absence of sur- 
face scattering. The constant of proportionality between 
5a j a and (A/£) 2 therefore lies between and — ir 2 in the 
long mean free path limit. Comparing this range to the 
values of \a\ stated above, we conclude that the magni- 
tude of the percent change in the effective conductivity 
due to roughness can be greater in either the long or short 
mean free path limits. 



V. GIANT MAGNETORESISTANCE 

Magnetic multilayers are alternating layers of ferro- 
magnetic and paramagnetic metals. For a large field in 
the plane of the layers, the magnetic moments of the fer- 
romagnetic layers align, creating the parallel magnetic 
configuration (P). At zero or a lower field, the magnetic 
moments of the layers are not aligned in parallel. With a 
proper choice for the thickness of the paramagnetic lay- 
ers, the ferromagnetic layers may be coupled antiferro- 
magnetically so that adjacent magnetic layers have mo- 
ments pointing in opposite directions. We call this the 
antiparallel (AP) configuration. The resistance of the 
multilayer in the two magnetic configurations is differ- 
ent, leading to a magnetoresistance. It is common to 
characterize this magnetoresistance as a ratio called the 
giant magnetoresistance or GMR, 



GMR = 



(37) 



where p AP and p p are the resistivities of the antiparallel 
and parallel configurations. 

In many systems the spin relaxation length, which is 
how far an electron's spin maintains its orientation, is 
long compared to the thickness of the layers. Hence, it is 
often possible to approximate the resistance as being due 
to two parallel channels for conduction, one channel for 
each of the two possible spin orientations. In the follow- 
ing we make this approximation, the corrections to which 
are well know in both the CPP and CIP geometriesc3E3. 



A. Long Mean Free Path Limit 

In order to deduce how the CPP GMR is effected by 
roughness in the long mean free path limit, within each 



spin channel we treat the layers as resistors in series. 
First, let tpM and tpM be the thicknesses of the para- 
magnetic and ferromagnetic layers, respectively, and let 
Ay = 2(tpM + tphi)- The resistivity for electrons of 
either spin in the paramagnetic material is denoted by 
Ppm, and the resistivity of the majority/minority elec- 



trons in the ferromagnetic material is Pfm,: 



laj/r, 



The 



interface resistance for the majority/minority electrons is 
P*maj/min Next ' define p ma j and p min by 



Ayp ma j 

Ayp 

l it i l l. 



2tpMPPM + 2tFMPFM (38) 
. Ill I II (39) 



In the parallel magnetic configuration, the resistivity of 
the majority spin channel is p ma j, and the resistivity of 
the minority spin channel is p m i n ■ Adding the resistivity 
of the two spin channels in parallel gives a resistivity for 
the parallel configuration of 



p PminPmaj 

Pcpp — ~~~ ■ 

Pmin i Prnaj 



(40) 



When the layers are aligned antiferromagnetically, both 
spin channels have the same resistivity, (p ma j + Pmin) I '2. 
Adding the two spin channels in parallel gives an antipar- 
allel configuration resistivity of 



Ap 1 / 

Pcpp — ~^\Pmaj 



Pmin ) • 



(41) 



As we have seen in the previous section, the effect of 
interface roughness in the long mean free path limit can 
be expressed as an increase in the interface resistance, 
5p*. To determine the effect of interface roughness on 
the parallel resistivity, the antiparallel resistivity, and the 
GMR, we expand Eqs. ©, ©, and @ to linear order 



in the changes in the interface resistances, 5p* ma ^ 
5p 

nun ' 



and 



{t FM + t PM )5p I 



2p 



2 

min 



maj ~r~ Pmin 
n maj 



2 ^ Pmaj 



iP: 

/ , \2 ^ fmin 

\Pmaj T~ pmin ) 



7$P*n 



(t FM + t PM )5p AP = -5p* maj + -5p* min 



(42) 
(43) 



5GMR = 



1 



1 



1 



Ay\ 

v Pmaj Pmin 



(44) 



5p 

min PminfiPmaj)- 



The changes in the interface resistances due to long 
length scale disorder may be obtained from Eq. ( |33| ) using 
the fact that p* is the inverse of a* for both the majority 
and the minority electrons, 



ma j j min 

o* 7, . 

maj j min 



5a 



maj fr 



maj fr 



(45) 



From Eq. ( |45D we can see that 5p* is positive for both 
the minority and majority electrons. Consequently, in 
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the long mean free path limit both Sp p and Sp AP are 
positive, i.e., they increase with increasing long length 
scale interface disorder. On the other hand, the second 
term on the right hand side in Eq. (jij) can be either pos- 
itive or negative depending on the values of the resistivi- 
ties and interface resistances. Hence when the electronic 
mean free paths are much greater than the layer thick- 
nesses, the CPP GMR can either be enhanced or reduced 
by roughness depending on the sample. The same holds 
true for the CIP GMR since in the long mean free path 
limit the CPP and CIP resistances are the same. 



B. Short mean free path limit 

To compute the GMR in the short mean free path limit 
we consider a four-layer repeat unit with layers 1 and 
3 being ferromagnetic and layers 2 and 4 being param- 
agnetic. As in the previous section, in the short mean 
free path limit we do not include interface scattering be- 
cause it does not change the results qualitatively. Let 
t fm = = Ay 3 be the thickness of the ferromagnetic 
layers and tpM — Aj/2 = Ay a be the thickness of the 
paramagnetic layers. The conductivity for either spin di- 
rection in the paramagnetic layers is defined as cfpm, and 
the conductivities for the majority and minority spins in 
the ferromagnetic material are <JpM,maj and (JFM,mim 
respectively. 

In the parallel magnetic configuration, the majority 
band electrons experience conductivities of o\ = 03 = 
&FM,maj and 1T2 = 04 = <JpM- Similarly, the minority 
band electrons experience conductivities of a\ = 03 = 
c™,mm and 02 = 04 = opM- When the magneti- 
zations of adjacent ferromagnetic layers are aligned an- 
tiparallcl, both spin channels have the same net conduc- 
tivity since the ferromagnetic layers (1 and 3) have the 
same thicknesses and the paramagnetic layers (2 and 4) 
have the same thicknesses. The conductivity for cither 
spin channel can thus be computed using ai = (JFM,maj, 
03 = OFM.mm, and <7 2 = (J4 = &PM ■ 

Solving Laplace's equation with the boundary condi- 
tions described in Sect. IV B, we compute numerically 
the conductivities for the parallel and antiparallel config- 
urations, op and (tap- The giant magnetoresistance then 
follows from Eq. (|37j), which when expressed in terms of 
the conductivity becomes {up— a ap) /tap- In contrast to 
the CPP conductivity in either the parallel or antiparal- 
lel magnetic configurations, which always increases with 
long length scale interfacial roughness, we find that the 
CPP GMR can either increase or decrease with (A/£) 2 
depending on the geometry of the sample and the layer 
conductivities. Specifically, for a sample with a given 
tpMi tpM, and <J f m ,maj / mim the GMR initially increases 
with roughness as Gpu is increased from zero. At some 
point below ((JFM,maj J raFM,min)/'2, <JpM reaches a criti- 
cal value at which the GMR is independent of roughness. 
The GMR subsequently decreases with roughness as upm 
is increased beyond this value. This critical value, <rp M , 



depends on the layer thicknesses and <JFM,maj/min and 
has the form 

°pm = C(l - ex P (-rt PM /Ay)), (46) 

where C is a constant slightly below (uFM.maj + 
CFM,mm)/2 and r is a constant of order 10, The exact 
value of these constants depends on <JFM,maj/min and the 
layer thicknesses. 

The CIP GMR in the short mean free path limit can be 
computed in the same manner as the CPP GMR by just 
changing the boundary condition from an applied field 
in the y-direction to an applied field in the x-direction. 
Here, we find the CIP GMR to be positive and propor- 
tional to (A/£) 2 . In the short mean free path limit, the 
CIP GMR therefore vanishes only when the interfaces are 
flat (A = 0). Figure | shows the CIP GMR as a function 
of crpM/vFM,min and a FM ,maj I^FM,min for a multilayer 
with A = 1.5, £ = 10, t FM = 5.33, and t PM = 2.67. The 
CIP GMR increases as tpm decreases, and within a par- 
ticular geometry the GMR reaches its maximum value 

when (TFM,maj ^S> UFM,min and (TpM ~ \&FM,maj- 




FIG. 9: Current-in-Plane GMR in the short mean free path 
limit. The presence of interfacial roughness leads to a giant 
magnetoresistance which is proportional to (A/t;) 2 and thus 
vanishes only for A — 0. The conductivities of the paramag- 
netic layer (for either spin), the ferromagnetic majority bands, 
and the ferromagnetic minority bands are upm, o-FM,maj, and 
c f M,min , respectively. The maximum GMR is achieved when 
tpM is small compared to tFM and for OFM.maj ~S> <JFM,min 
and apM ~ ^<JFM,maj- The data shown correspond to a mul- 
tilayer with A = 1.5, £ = 10, tFM = 5.33, and tpM = 2.67. 

In addition to solving for the conductivities exactly 
in the short mean free path limit, one can also use the 
approximate expressions of Eqs. (p6|), (|29|), and (|32"|). Al- 
though not as accurate as the numerical solution, this 
technique may be useful in estimating the size of the ef- 
fect for an actual experiment. 

According to these equations the percent change in the 
conductivity due to long length scale surface roughness 
is proportional to the ratio of the CPP and CIP conduc- 
tivities corresponding to flat interfaces. There are two 
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possible magnetic orientations of the layers, parallel and 
antiparallcl. For each of these magnetic orientations the 
electrons can either be in the minority or majority bands. 
Thus, since there are two geometries (CPP or CIP), two 
magnetic orientations (P or AP), and two spin channels 
(maj or min), there is a total of eight conductivities to 
be specified. We denote the conductivity of the majority 
band for the CPP geometry in the parallel magnetic ori- 
entation as crQpp ma j and label the other conductivities 
in a similar fashion. 

For the CPP geometry the conductivities can be com- 
puted by adding the resistances of the layers in series. 
Using the notation for p ma j and p m i n in Eqs. ( |38| ) and 
(E39l) with p* = 0, the CPP conductivities are 



below. 



p 

fJ CPP,maj 

p 

a CPP,min 

AP _ AP 

a CPP,maj — a CPP,min 



1/ Pmaj 

i/a 



mm 

2 



Pmin Pmaj 



(47) 
(48) 

(49) 



For the CIP geometry the conductivities can be com- 
puted by adding the resistances of the layers in parallel. 
It is useful to define the analog of Eqs. ( |38| ) and ( |39| ) for 
the layers in the parallel case, 



Aycr maj 
Aya 

mtn 



= 2b 



1 



PM- 



2b 



1 



FM " 



PPM PFM,maj 
1 1 

2bpM V 21fm ■ 



PPM 



PFM.i 



(50) 
(51) 



Note that with this notation p ma j l/a ma j, which will 
be important below. The CIP conductivities are then 



Spcpp 



SPcpp 



$PCIP 



X AP 
OPCIP 



P 2 ■ 

1 m 1 n 



( a maj Pmaj) 



(Pmin H~ Prnaj)^" 



(55) 



j \P min 



{Pr, 



~t~ Pmaj)^ 



= 16 - 



16 - 



= 16 - 



&CIP 



1 



(o-c/p) 5 



(o"c/p) 2 



Pcpp 



&CIP 



&CIP 



(56) 

-^—\ (57) 
Pcpp > 



1 



Pcpp 



(58) 



Using the definitions of p ma j, p m in, &7na 3 , and a mtn , it 
follows that in the CPP geometry the change in resistiv- 
ity is negative, Sp ma j/ min < 0, while in the CIP geometry 
the change in the resistivity is positive, $Pmaj/min > 0. 
This result is expected from Figs. || and |[ In the CPP 
geometry the waviness of the interfaces allows the cur- 
rent to travel a greater distance through the less resistive 
layers, thereby reducing the resistance. In the CIP ge- 
ometry the waviness of the interfaces disrupts the current 
flow through the low resistivity layers, increasing the re- 
sistance. 

The GMR is determined by the ratio of the parallel 
and antiparallel resistivities. As seen above, the paral- 
lel and antiparallel resistivities either both increase or 
decrease depending on the geometry. Thus, one must 
compute their ratio explicitly to determine whether the 
GMR increases or decreases with long length scale inter- 
face disorder. Expanding Eq. ( |37j ) for the GMR to linear 
order in Spp and Spap and substituting the changes in 
the resistivities of Eqs. (55) - (|58|), the change in the 
GMR for the CPP and CIP geometries is 



a CIP,maj 
P 

"CIP,min 

AP _ AP 

a CIP,maj — a CIP,min 



&CIP 



(52) 
(53) 

(54) 



Within the two channel conduction model, for a given 
geometry and magnetic orientation the conductivity of 
the sample is the sum of the conductivities of the two 
spin channels. According to Eqs. (52) - (|54|), within this 
model the parallel and antiparallel configurations in the 
CIP geometry have the same sample conductivity, acip- 
Due to long length scale interface disorder each of the 
conductivities in Eqs. (E^) - (|54|) changes according to 
Eqs. ((2^), (p9|), and (|32|). The net conductivity and re- 
sistivity of the sample therefore also change. As in the 
previous section we denote the change in the resistivity 
by Sp — —Sa/cr 2 . Here, the resistivities are different for 
the two geometries as well as the two magnetic orienta- 
tions, so there are a total of four <5p's, which are given 



SGMR 



CPP 



AP 

Pcpp 



SGMRcip 




Using Eqs. ( |59| ) and (|60|), it can be shown that in the 
short mean free path limit the GMR decreases in the CPP 
geometry (SGMR C pp < 0) and increases in the CIP ge- 
ometry (SGMRcip > 0). In the CIP case, this result is 
consistent with our numerical computation of the GMR, 
which always increased with long length scale interface 
roughness. In the CPP case, however, we found numer- 
ically that the GMR could either increase or decrease 
with interface roughness depending on the sample geom- 
etry and the conductivities of the paramagnetic and fer- 
romagnetic layers. In both geometries, Eqs. (^9j) and ( |60| ) 
generally yield much larger changes in the GMR than we 
predict numerically. This discrepancy is due to the fact 
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that the agreement between the analytic expressions for 
a and the numerical values varies depending on the layer 
conductivities. If the analytic and numeric values for a 
agree very well within a particular spin channel, then the 
analytic expression for the change in the conductivity for 
that channel will agree very well with the change com- 
puted numerically. When one considers a different spin 
channel with different layer conductivities, the analytic 
a will generally be a worse approximation of the numeric 
value because the layer conductivities have changed. The 
change in the conductivity for this channel computed an- 
alytically will likewise be a worse approximation of the 
numeric value. This will result in deviations between the 
antiparallcl and parallel conductivities computed analyt- 
ically and numerically. Since the GMR involves a ratio of 
these quantities, the error will be more significant in the 
GMR, resulting in large deviations between the analytic 
and numeric changes in the GMR. 



C. Estimates 

Both the long and short mean free path effects de- 
scribed in this paper must be present to some extent in 
any sample. The real question is how large these effects 
are and whether they can account for what is seen ex- 
perimentally. In this section we estimate the size of the 
two effects in Fc(3nm)/Cr(1.2nm) multilayers using ex- 
perimentally determined parameters. Although there are 
no adjustable parameters in these calculations, they are 
still merely estimates because the actual experiments are 
in neither the long nor the short mean free path limits, 
but somewhere in hetween. This can easily be seen from 
Cyrille et aZ.'s datacd. In the long mean free path limit 
the CPP and CIP resistances are the same, whereas ex- 
perimentally the ratio of the CPP to the CIP resistances 
is roughly 1.5. In the short mean free path limit there 
is no CIP GMR for flat interfaces, while there is roughly 
a 10% CIP GMR observed experimentally. In addition 
to taking the long and short mean free path limits, the 
calculations in this paper are performed in two dimen- 
sions instead of three dimensions, and we use one of the 
simpler Boltzmann equations. 

To estimate the size of the long mean free path effect, 
we use Eqs. ( ^2| ) and (|43|). The resistivities, p m in and 
Pmaj, may be determined experimentally from Cyrille et 
al.'s measurements of p p and p AP using Eqs. @ and 
( |il| ) . Although the data vary with the number of bilay- 
ers, the results depend only weakly on which data points 
one uses to compute p m aj and p m in , with the values in the 
range p m in — 45.7 ±0.9 pflcm and p m aj = 130±8/iS!cm. 
These are 'very close to the values obtained from Zam- 
bano's datafij of p m in = 45/if2cTO and p m aj — 143/iQcm, 
which is not surprising because the overall resistances 
in the two sets of experiments are close together, even 
though the trends are not the same. 

To get the changes in the majority and minority sur- 
face resistances, 5p* ma j and 5p* min , one needs to know 



the values of the surface resistances that the changes will 
be computed from. These parameters have not been de- 
termined in the experiments by Cyrille et al; however, 
they have been determined in a series of experiments by 
Zambano et al. in which measurements were taken, with 
different layer thicknesses and number of layersEil. In 
their notation p* ma j and p* min are equal to AW F& , Cr and 

ARp e ^ Cr , respectively, and they find these values to be 

Pmaj = 2.7/fhn 2 and p* min = 0.5/Slm 2 . We will use 
these numbers for our estimate because the overall re- 
sistances in the two sets of experiments are close to one 
another. 

Finally, we need to determine the waviness or rough- 
ness of the interfaces in the experiments. This was 
quantified by Cyrille et al. using two techniques - low 
angle x-ray diffraction and energy filtered imaging us- 
ing cross-seetioned samples in a transmission electron 
microscope.Eil Both techniques show that the root mean 
square deviation of the layer height increases within the 
multilayer and has the form a = a a M a , where a a — 
0.37nm is the root mean square deviation of the first bi- 
layer, M is the bilayer index, and the exponent a equals 
0.4. For the samples deposited at a constant pressure of 
5mTorr, the average distance between surface bumps, or 
the period £ in our model, was determined to be lOnm, 
independent of the bilayer indexfij For a multilayer with 
N bilayers, the average value of (Aft;) 2 used in Eq. ( |45| ) 
to compute the changes in the interface resistances is thus 



1 N 



(61) 



M=l 



where we have used the fact that the root mean square 
fluctuation for a sine wave is A/y2. 

The result of using Eqs. @, (g§), @, and @ to es- 
timate the size of the long mean free path effect is shown 
in Fig. [l0| As expected both the parallel resistivity (solid 
line) and the antiparallel resistivity (dotted line) increase 
with roughness. The size of the increase in the antiparal- 
lel resistivity measured in the experiments by Cyrille et 
al. is comparable to this estimate for the antiparallel re- 
sistivity. However, our estimate of the size of the increase 
in the parallel resistivity is much larger than the negligi- 
ble or possibly negative change in the parallel resistivity 
seen in the same experiments. 

In the short mean free path limit, there are four re- 
sistivities for the two geometries (CIP and CPP) and 
magnetic configurations (P and AP). The changes in the 
resistivities for these four cases are given in Eqs. (|55| ) 
- (58). In evaluating these equations, we use the CPP 
resistivities, p p PP and p^pp, determined by the same 
Pmaj and Pmin used in the long mean free path limit us- 
ing Eqs. @ and @. The CIP conductivity , a C ip, is 
taken to be in between the observed parallel and antipar- 
allel conductivities: {&cip)~ 1 — 25pflcm. Finally, the 
conductivities u ma j and <T m i n are not readily determined 
from experiment. To estimate them, we impose the con- 
dition a maj /a min = Pmm/Pmaj and use the value for 
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path limits. It is entirely possible that for intermediate 
mean free paths the two effects cancel, leading to the neg- 
ligible change in the resistivity seen in the experiments 
by Zambano et al. If the long mean free path effect does 
dominate, then a result similar to that of Cyrillc et al. 
would be observed. In the CIP configuration there are 
no complications due to the effects canceling. The esti- 
mated change in the resistivity in both limits is positive 
and comparable to the increase observed by Cyrille et al. 



N 



VI. CONCLUSIONS 



FIG. 10: Estimates in the long and short mean free path lim- 
its for the changes in the parallel and antiparallel resistivities 
due to long length scale interface roughness. The expressions 
for the changes in the resistivities are evaluated using experi- 
mentally determined parameters, but the data shown are only 
estimates because the experiments are in between the long 
and short mean free path limits. The magnitudes of the resis- 
tivity changes are consistent with those measured by Cyrille 
et al; however, for intermediate mean free paths there can 
be cancellations of Sp p ^ AP and &p^Jpp ', which may explain 
the results of Zambano et al. For the CIP geometry there is 
no such cancellation and the estimates are comparable to the 
changes seen by Cyrille et al. 



a cip = &min + &maj from above. The resulting changes 
in the resistivities are shown in Fig. [l0| As expected, 
here the changes in the resistivities for the CPP config- 
uration are negative, while the changes in the resistivi- 
ties for the CIP configuration are positive. Although it 
was mentioned previously that the analytic expressions 
for the changes in the GMR given in Eqs. ([59|) and ( J60| ) 
predict a much larger effect in the short mean free path 
limit than we predicted numerically, the changes in the 
resistivities presented here are more reliable. Since the 
GMR involves a ratio of the antiparallel and parallel re- 
sistivities, the error in the analytic GMR will be more 
significant than the error in the analytic antiparallel and 
parallel resistivities. 

As Fig. |l0| illustrates, the changes in the CPP resistivi- 
ties are of roughly the same magnitude, but opposite sign, 
in the long and short mean free path estimates. Thus, 
in a sample with intermediate mean free paths, there 
may be substantial cancellation of these two effects. For 
the CIP configuration, however, both the long and short 
mean free path effects tend to increase the resistivity by 
a similar magnitude. An increase in the CIP resistivity 
of comparable size to these estimates is indeed seen in 
the experiments by Cyrille et al. 

These estimates show that the effects described in this 
paper are of the correct size to describe the observations 
of Cyrille et al; however, detailed quantitative compar- 
ison is not possible because the mean free paths in the 
experiments are in between the long and short mean free 



In this paper we have examined the effect of interface 
disorder which is long on the atomic scale. These kinds 
of fluctuations are ubiquitous in metallic multilayers. A 
semiclassical Boltzmann equation was solved in both the 
limits where the electronic mean free paths were short 
and long compared to the layer thicknesses. 

In the short mean free path case the current flow is 
nonuniform, and long length scale interface disorder in- 
creases the effective conductivity in the CPP geometry 
and decreases the effective conductivity in the CIP ge- 
ometry. In the CPP case, the effect is due to interface 
disorder providing a less-resistive, non-linear path of cur- 
rent flow. In the CIP case, the effect results from inter- 
face roughness disrupting the flow of current through low 
resistance layers. In the long mean free path case the 
current flow is uniform in the CPP and CIP geometries. 
The resistance increases with long length scale interface 
roughness in both geometries because of the additional 
scattering created by longer interfaces for the disordered 
layers than the flat layers. 

The experiments discussed in this paper are in nei- 
ther the short nor the long mean free path limits, but 
somewhere in between. Nonetheless, in estimating the 
size of the long and short mean free path effects we find 
that in the CPP geometry the increase in the antipar- 
allel resistivity observed by Cyrille et al. is comparable 
to the increase we predict in the long mean free path 
limit. The increase estimated in the CPP geometry for 
the parallel resistivity, however, is larger than the neg- 
ligible increase measured by Cyrille et al. We find that 
the long and short mean free path effects tend to cancel 
for intermediate mean free paths in the CPP geometry, 
which could explain the observations of Zambano et al. 
For the CIP geometry, in both the long and short mean 
free path limits we estimate an increase in the resistivity 
that is comparable in magnitude to the increase seen by 
Cyrille et al. Therefore, the effects described here may 
be the source of the experimentally observed increase in 
the resistivity with long length scale disorder, although 
clearly more theoretical investigation is needed to under- 
stand the crossover between the long and short mean free 
path limits. 
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